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Abstract 

In hep-th/0501082, a field theoretic "toy model" for the Landscape was proposed. We 
show that the considerations of that paper carry through to realistic effective Lagrangians, 
such as those that emerge out of string theory. Extracting the physics of the large number 
of metastable vacua that ensue requires somewhat more sophisticated algebro-geometric 
techniques, which we review. 



1. Introduction and Summary 



One of the striking observations of KKLT P was that the existence of a large number 
of long-hved metastable vacua in certain compactifications of string theory (dubbed the 
"landscape" [21) provides a concrete instance in which the anthropic principle |3|4|5|6|7|8|9j 
might be realized in Nature. Unfortunately, many discussions of these ideas get bogged down 
because it is hard to disentangle the intricacies of the string theoretic constructions from the 
"anthropic" questions one would like to address. 

The paper of Arkani-Hamed, Dimopoulos and Kachru ^Oj was, therefore, very useful 
in clearing away the string-theoretic underbrush and presenting a simple, tractable field- 
theoretic model with a large number of vacua in which anthropic questions could be addressed 
(see also [H]). 

One of the concepts to emerge from their investigation was the notion of a "friendly 
landscape." In general, all of the couplings, Ca of the theory will vary between the different 
vacua. However, there is a qualitative difference between those couplings which "scan" (those 
whose standard deviation is much larger than their mean value) and those couplings which 
"don't scan" (those which are sharply-peaked about their mean value). In making anthropic 
arguments, one usually considers the situation in which one coupling is allowed to vary, while 
the others are held fixed. If all the couplings vary appreciably, the anthropic bounds are 
much weaker, or go away entirely. So it was very useful for the authors of jTOj to identify a 
mechanism by which the couplings one would like to "tune" anthropically scan, whereas the 
remaining couplings are sharply peaked. 

Their model had three basic features: 

1) A large number, A^, of scalar fields, 0,. 

2) A decoupled form for the scalar potential, V{(l)) = J2i Vi{4>i)- 

3) A decoupled form for the 0-dependence of the observable couplings of the model. 

The first has a fairly natural realization in string theory. Many compactification of string 
theory have hundreds of moduli which, when the physics which lifts the vacuum degeneracy 
is included |12li;-ill4ll5lieill7ll8ll9l2()l21l22j . form natural candidates for the 0^. We can, 
quite plausibly, take "A^" to be the number of (complex) moduli. However, 2,3) are rather 
unnatural from this point of view and, indeed, it's hard to imagine that such a decoupled 
structure might emerge from string theory. 

In the present work, we would like to overcome this drawback and present what we hope 
is a realistic version of the scenario of ^U]. Our approach will be to start with, essentially, the 
most general A/" = 1 supersymmetric effective field theory with a large number of "moduli" 
chiral multiplets, $*. We will then see what conditions must be imposed in order to realize 
a friendly landscape. 

In ^ we will discuss the theory of A^ chiral multiplets coupled to A/" = 1 supergravity. 
Whereas any such theory must be cut off at a scale M^. < Mp, we will see that, at large- 
A^, we will will need to impose the stronger condition Mc < Mp/ \/N, in order to have a 
sensible effective field theory. Moreover, we will find, in ^2.H that the couplings among the 
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moduli chiral multiplets must be suitably small; we will summarize the condition that we 
must impose on these couplings by saying that they are generically small. We then turn to 
a polynomial truncation of the superpotential of the model. While not indispensable, such a 
truncation makes the analysis of the vacuum structure amenable to perturbation theory. As 
usual, for this to be a reliable guide to the vacuum structure, the super-renormalizable terms 
in the superpotential must have coefficients governed by a mass scale M^, <C M^. f ^2.2p . In 
mi.'iV we will see how these considerations mesh with the most popular arena for landscape 
considerations — F-theory vacua with fluxes. 

Next ( §2.4j) . we impose the discrete R-symmetry found by ^01 to lead to a friendly 
landscape for the cosmological constant, and use the GL{N, C) symmetry of fleld redeflnitions 
to simplify our problem f ^2.5p . 

In ^ we lay out the algebro-geometric techniques used to determine the vacuum structure 
and extract information about the distribution of values for various physical quantities among 
the 2^ vacua of the theory. In ^ we use these techniques to compute certain "holomorphic 
moments" which characterize the distribution of values of the cosmological constant. In 
S|5J we discuss the physics that ensues when one assumes that the couplings are chosen from 
some (unspecifled) probability distribution, generahzing the considerations of [TU] . We apply 
our analysis both to the superpotential and to other holomorphic couplings. 

Finally, in ^Hl we discuss some generalizations of our techniques and future directions. 



2. General Features of SUSY Landscape Sectors 

As a fleld theoretic model for the landscape sector, we could start with an arbitrary SUSY 
fleld theory of chiral superflelds and A^' vector superflelds, all coupled to A/" = 1 SUGRA. 
For simplicity, we will restrict ourselves to the case A^' = and only briefly touch on gener- 
alizations to gauged hidden sectors in this work. Thus, the vacuum structure of the model 
can be determined using the two-derivative effective action of the SUSY non-linear cx-model 
describing the A^ chiral flelds at energies below a cutoff scale, Mg. 

Clearly, for any sort of effective fleld theory to be valid, we must take Mc < Mp. However, 
as noted in ^U], when one has a large number, A^ ^ 1, of flelds, radiative stability of Newton's 
constant requires 

1 

m| < N P^^' 

To see this, note that the action for A/" = 1 supergravity interacting with A^ chiral flelds is, 
in superspace notation (we use Mp for the reduced Planck mass and the conventions of j23j). 

S = -3M2 j d^zE-^e'^^^'^'^^ + j (fz(t)^W{^) + h.c. (2.2) 

where is the superdeterminant of the vielbein, 2; is a superspace coordinate, and is a 
compensator superfleld. The Einstein-Hilbert term comes from the leading ($- independent) 
piece of the flrst term in (j2.2p . At one- loop, this receives a quadratically-divergent contribu- 
tion, 



$=$=0 



(i§)w-iv(3§) (2.3) 
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where g^j = didjK is the Kahler metric of the cr-modeL The enhancement which comes from 
having fields running around the loop requires us to set the cutoff, to be parametrically 
smaller than Mp. 

In A/" = 1 supergravity, the chiral multiplets parameterize a Kahler manifold, whose 
Kahler form, uj = ^ddK. The superpotential, ly ($) transforms as a section of a line bundle, 
C, whose first Chern class, 

cJC) = l-u = — ^—^ddK. (2.4) 

The fiber metric on £ is 

h{^, $) = e^(*'*)/*^^' (2.5) 

and the connection is of type (1,0), Di = di + diK/Mp. Under Kahler transformations, 
^($) ^ e^(*)/*^pV(<l>), 

ir($\ ¥) K{¥, ¥) + /(<l>^) + f{¥), (2.6) 

The supersymmetric vacua of this model are the critical points of the superpotential with 
respect to the Chern connection, i.e. the points at which DiW e T(T*X ® C) intersect the 
zero section of T*X ® C, 

DiW = diW + {diK)W/Ml = 0. (2.7) 

More generally, all the vacua of this model, supersymmetric or not, are critical points of its 
scalar potential, 

V = e^l^^l [g'W.WDjW - SII^IVM^) . (2.8) 

Since we will be interested in vacua in which 1 1^ | ~ < -C , we will always be safe 
in neglecting the connection term in ()2.7p and hunting for ordinary critical points, diW = 0. 

In fact, to justify a perturbative analysis of the effective Lagrangian (j2.2j) . we require 
that H^($) have a convergent Taylor expansion in a polydisk, |$j| < Mc- Just as the 
radiative stability of Newton's constant imposed constraints on the cutoff scale, Mc of our 
effective Lagrangian, radiative corrections to the Kahler metric imposes, at large-N, further 
constraints on the coefficients. 



2.1. Quantum Corrections to the Kahler Potential and Large Scaling 

Before embarking upon a study of the vacua of these models, we should study the radiative 
stability of these models in the limit of large A^. As our theory is an effective field theory, 
we expect that the characteristic radius of curvature for X is given by the cutoff scale, 
Mc- Thus, in a small enough neighborhood of a smooth point p G X, we can choose local 
coordinates such that $*(p) = and the Kahler potential takes the form (modulo Kahler 
transformations) , 
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where the Ki^i_ are dimensionless and symmetric in the multi-indices I„ and In and Qij is 
the Kahler metric at p. Further, as the superpotential W is holomorphic, it can be locally 
expanded as a power series in the holomorphic coordinates about the origin = 0, 

= Ao + y^ . . . = M^Wo + V , . . . (2.10) 

n,In 

where the are symmetric in the multi-indices I„ and have mass dimension 3 — n while 
the = MJ}~^Ai^ are dimensionless. 

While the superpotential is holomorphic and not renormalized, the Kahler potential is 
renormalized. It is important to check the radiative stability of its assumed form given the 
form of W in the large N limit. More precisely, we compute the effective Kahler potential 
(in the Wilsonian sense) at a scale Mc' lower than Mc by integrating out the modes in a shell 
of momenta between Mc and Mc', and require that the corrections are not parametrically 
larger in than the bare values. As we will see, this requirement will restrict the asymptotic 
growth of both the Aj^ and Ki^i_ parametrically in A^. 

Let us first consider the one-loop corrections to the Kahler potential coming from the 
superpotential W, 
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(2.11) 



where denotes the uncorrected inverse Kahler metric at the point p, the origin in field 
space. For n = n = 1, we get a one loop correction to the Kahler metric at the origin p 
arising from the dimensionless superpotential couplings Aij^, 

- log ^) 9''9'%MA,,r. (2.12) 

Radiative stability requires that Agfj^'' is "small" compared to gij, that the coordinate- 
invariant norms of tangent vectors to X at p under A^fj-^-* are smaller that those under 
Qij, or roughly, 

g^^Agff ~ log g^^g'^YAkiAf^i < g^^g.j = 5^ (2.13) 

We can re-phrase this requirement in precise, coordinate invariant terms by taking traces 
and determinants of both sides, 



det 



g^^Agff <1, Tr^^^A^^.V'^ < A^. (2.14) 



As promised, this gives us parametric bounds on the growth of the dimensionless superpo- 
tential couplings Aijk, 

g''g"g''A,,Aj-, ~ 0{N), (2.15) 



det ij 



9 9 9 ^zki^jki 



0(1). (2.16) 
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Since Aijk is a tensor, the interpretation of these bounds on the size of its components is, of 
course, a coordinate dependent question. Now, in a Kahler manifold, du = imphes that it 
is always possible (see |23] p. 107) to choose local holomorphic coordinates about any smooth 
point p such that = and 

didjK{^, <l) = Qiji^, $) = 5ij + terms of order > 2 in the $, <l. (2.17) 

In these coordinates the above bounds simplify to. 



Y,\A^Jk\^ -OiN), det 



jkl 

. k,l 



0{l). (2.18) 



Thus, in these special coordinates, we see that in order to satisfy these bounds, the generic 
components of Aijk must have parametrically suppressed magnitudes at large iV, 

\A,,k\ ~ 0{N-^) (2.19) 

though 0{N) of them can still be as large as 0{1) in that limit. Therefore, we will refer to 
a tensor Aijk obeying ()2.15|1 as being generically small. 

We can similarly obtain bounds on the coefficients Ai^ and Ki^i^ by considering the 
leading loop correction to the quadratic term in the Kahler potential coming from the cor- 
responding higher order terms. One finds a slew of conditions similar in form to ()2.18p . For 
instance, 

^2n-2^nn . . ^^^^^^^ ^^^^ ^ o{N), (2.20a) 

det [Mf -2^7^ V^'^ ■ ■■g'-^''-^'A,...,„^,iAi,..,„^,j\ ~ 0(1), (2.20b) 

g^^'^...g^'^'-K,,...,„,,...,^^0{N), (2.20c) 
det [g^^g'^'^ ■ ■ ■ ^ 0{1). (2.20d) 

In particular, this implies that 0{N) of the coefficients Ai^ could be 0{1), while the generic 
coefficient must be small, 

|Ai...iJ ~ M-('^-=^)0(iV^(""^)/2). (n>3) (2.21) 

Thus, just as with Aijk, we will refer to any tensor obeying bounds of the form ()2.20|) as 
generically small. In general, these conditions, as well as a very large number of others 
corresponding to loop diagrams involving multiple vertices, constrain all the coefficients of 
the higher order terms in the Kahler potential and superpotential to be generically small at 
large A^.^ 



^ We note that there is an additional constraint on the size of the coefScients coming from the requirement 
that the power series we have been writing down actuaUy absolutely converge in a region of size ^ Mc- These 
constraints, while stronger than those coming from radiative stability, only constrain the asymptotics of the 
Ai^ for large n, rather than constraining the terms at some particular order in n. 
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2.2. SUSY Vacua of the Renormalizable Wess-Zumino Model at Large 



We now wish to study the vacuum structure of this effective field theory. As discussed in 
the previous subsection, this amounts to studying the critical points of the superpotential, 
diW = 0. For large N, we expect that the number of critical points of W within a polydisk 
of radius Mc to scale exponentially with A^. 

Since we assume that the power series expansion for W is absolutely convergent in the 
polydisk, a field theorist might reasonably take the approach of truncating the power series 
at some finite order and looking for the critical points of the resulting polynomial. Of course, 
we cannot hope that a truncation to any finite order polynomial can be an accurate guide 
to the critical points in the entire polydisk of radius M^. At best, we will hope to obtain the 
critical points inside some much smaller polydisk, of radius -C M^. Generically, however, 
we would not expect to find any critical points inside this smaller polydisk. The criterion for 
finding (trustworthy) critical points within this smaller polydisk for small A^ is well-known: 
we require that the coefficients of the linear and quadratic terms in ()2.1()|1 be Ai-^ ~ 0{M!^) 
and Ai^ ~ 0{Mr), respectively, where 

Mr 

e < 1 (2.22) 



M, 



c 



As long as the quartic and higher terms in the superpotential are not anomalously large, 
they give negligible corrections to the critical points determined by truncating W to cubic 
order^. 

At large A^, this is not quite sufficient. Even assuming that the higher Aj^ are generically 
small, in the sense of the previous subsection, we still need to require 

e < ^ (2.23) 



in order for these higher-order corrections to be negligible. To see this, we can look at the 
invariant quantity, 

iM;^g'-^d,WdjW)\^=^^, (2.24) 

where is the critical point derived from the cubic approximation of W. Using the generic 
smallness of the \AiJ ~ Mc"^""^^Ar-("-i)/2 and l^^,] ~ M, = eM^, we see that the corrections 
grow parametrically with A^, unless ()2.23|) is satisfied. In what follows, we will keep e as a 
free parameter, cognizant of the fact that it must be sufficiently small for the story to work. 

In non-supersymmetric theories, ensuring that under radiative corrections the coefficients 
of the super-renormalizable terms in the potential remain much smaller than the cutoff 
is called the hierarchy problem. In a supersymmetric theory, there are no perturbative 
corrections to the superpotential. Nonetheless, the fact that the coefficients of the super- 
renormalizable terms in the superpotential are of order Mr, rather than Mc, means that 
the point about which we are expanding is, in some sense, "special." That will be further 
brought home in §2.4^ where we will assume that the point $ = will be a point where there 
is an unbroken Z4 R-symmetry. 

^As long as we are away from the discriminant locus to be discussed below 
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2.3. Distributions of Models and Flux Compactifications of IIB 

Since F-theory compactifications witli flux are tlie prime motivating example of "landscape" 
models with a large number of vacua, let us pause to consider how such models fit in with our 
general considerations, as developed so far. Any given set of fluxes satisfying the requisite 
tadpole cancellation conditions in an F-theory or IIB Orientifold compactiflcation gives rise 
to a Gukov-Vafa-Witten [23] superpotential for the complex structure moduli. Since there 
are, in general, many solutions to the tadpole cancellation conditions, we have an ensemble 
of theories with different superpotentials, labeled by the possible fluxes. 

In the IIB case, the elements of the ensemble are labeled by all choices of integer fluxes 
through the 26^'^ + 2 three cycles of the Calabi-Yau 3- fold Ai, 

^ ^ F = iv, e Z, —\— f H = MaeZ. (2.25) 



(27r)Vys„ (2vr)V 



compatible with the tadpole cancellation condition for induced D3-brane charge on Ai, 



^3 = ttJt^ [ FAH. (2.26) 

'M 



(27r)4a 



where includes the contribution of space-fllling, mobile D3-branes, D7-branes, and 03- 
planes. Now, if Q3 ~ 62,1 ~ 0{N) then tadpole cancellation requires that each of the integer 
fluxes can then be at most Na ~ Ma ~ 0{N^^'^). Each choice of fluxes gives rise to a Gukov- 
Vafa-Witten superpotential for the IIB axio-dilaton r and the complex structure moduli 
Zi = J^^^oi M, 

W = [ GAn. (2.27) 
Jm 

where G = F — tH and the Ai and form a symplectic basis of 3-cycles {Sa} = {Ai, B^} 
in A4 with respect to the intersection pairing on H^^M.^'L). In particular, we see that the 
tadpole cancellation condition places bounds on the growth of W on N , 

WdJ^ [ G I V(iV^ - TM')zi ~ 0{N). (2.28) 

which is certainly consistent with the condition that the couplings of the moduli are indeed 
generically small. However, it is not at all clear that there is any point in the moduli space 
about which the cubic truncation gives a good approximation to the locations of its critical 
points - that is, we do not expect the renormalizable approximation to be a useful guide to 
the vacuum structure. 

Further, notice that the superpotential is odd under G —G, while the tadpole can- 
cellation condition is even, so if G satisfles the tadpole constraint, then so does —G. This 
means that if we flnd a supersymmetric vacuum z*, a critical point of ()2.27|) with a given 
flux G, and a corresponding value W* of the superpotential, then z* also corresponds to a 
vacuum of the model with flux —G and superpotential — W*. Thus, even though the average 
value of W at the supersymmetric vacua for any fixed choice of the fluxes may not vanish, 
the ensemble average of W certainly does vanish. 
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Motivated by this example, we will be interested in situations where we are actually given 
a distribution or ensemble of models. That is, we will assume that physics at high energies 
~ Mp can be understood as providing a distribution of coefficients Aj^ for the effective 
low-energy < Mc models described above. Indeed, in the IIB flux vacua, the physics which 
determines the values of the quantized fluxes is not even field theoretic in nature - it likely 
involves high energy string/brane dynamics, topology change, etc. In particular, the scanning 
of the cosmological constant in these models is explained by high energy physics. Certainly, 
we could accept such a high energy explanation and restrict our consideration to distributions 
of coefficients which are symmetric under Ai^ — >■ —Ai^ and share this property of flux vacua. 
However, we will focus on situations where there may be a low energy explanation for this 
scanning.^ One way to achieve this ^01 is through the imposition of an R-symmetry, which 
we discuss presently. 

2.4. R-Symmetry, The Renormalizable Wess-Zumino Model, and A 

The general renormalizable Wess-Zumino model of interacting chiral superfields is 
described by a quadratic Kahler potential, 

K{<!>,<^)=g,jm\ (2.29) 

and a cubic superpotential, 

W{¥) =A + Ai(l>' + -Aij$'<l>^' + (2.30) 

2 3 

Radiative stability of the scaling of the Kahler potential at large restrict us to the case 
that generically \Aijk\ ~ A^~^ with 0{N) terms which are 0{1). The supersymmetric vacua 
for this theory are points where, 

diW = A + A,^^ + Ai,fe<l>^$^ = 0, (2.31) 

a set of A^ complex algebraic equations in the $*. In particular, these polynomials determine 
an ideal^ in the ring of polynomials in the (diW) C C[$^, . . . , $^]. Bezout's theorem 
(see Chapter 3, Theorem 5.5 of [20]) guarantees that for a generic choice of the coefficients 

•^Somc high energy input, however, may be inevitable. For example, one might worry about the origin of 
the small parameter e = M^/Mc that we required in order to make the cubic approximation. This appears 
to be a tuning of 0{N) relevant couplings in the model. Without some further explanation, one might worry 
that, together, these represents a fine-tuning of order e^, which would wipe out whatever "advantage" we 
gained in having 2^ vacua. This need not be the case if the smallness of each of the 0{N) couplings 
has a common explanation. Perhaps there is an approximate symmetry, broken weakly by the effects that 
generate the superpotential, which guarantees that these couplings are small. Indeed, such a symmetry 
would be reflected in symmetries of the resulting probability distribution for the couplings. For example, if 
the coefficients of the relevant operators were all selected from the same distribution, we would indeed only 
require a single fine tuning of the distribution. 

^See section O for a basic review of the commutative algebra language used here, and |26I27| for a more 
detailed introduction. 
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of W,^ these simultaneous quadratic equations have 2^ roots. From an algebraic point of 
view, this translates into the fact that the quotient ring, 

C[$\...,$^]/(9,iy) = C^^ (2.32) 

is a 2^-dimensional vector space over C (see Chapter 3, Theorem 6.2 of [2^1) generated by 
the images of the monomials ■ ■ ■ ^^'^ with ii < ■ ■ ■ < ij. < N. Explicitly, the quotient ring 
should be understood as the "polynomial functions" on the algebraic variety cut out by the 
equations generating the ideal, which in this case are the functions on a set of 2^ points, i.e. 
a 2^ dimensional vector space. 

Further, we saw that the critical points are generically contained in a neighborhood U 
of radius ~ M,. about the origin. Note that if we assume that the constant term A in is 
also of order < NM^, the superpotential at each of the critical points in U is roughly of the 
order of^ 

W ^ X 0{N). (2.33) 
This gives a supersymmetric contribution to the vacuum energy 

A = -3^ (1 + 0{M^/M^)) ~ -Ne^M^, (2.34) 
p 

where e = Mr/Mc < and M^/Mp < 1/^/N. Note that all corrections coming from 

the inclusion of the Chern Connection terms and higher order terms in the Kahler potential 
and superpotential are parametrically suppressed at large A^. 

When supersymmetry is broken, we get a positive contribution, A5 to the vacuum energy. 
The hope is that the distribution of values for W among the 2^ vacua will be such that (j2.34j) 
can very nearly cancel Ag.'^ To achieve this ITU], one wants a distribution of values of W , 
such that the standard deviation is large, compared to the mean value, ((VT)). 

As the authors of jTUj noted, this is easy to arrange. If the superpotential is odd under 
— then supersymmetric vacua come in pairs $** and — $** with opposite values 
of W and we would have ((W^)) = 0. We can enforce this by imposing a Z4 R-symmetry 
$*(?/, 6*) —^^{y,i6), under which the superpotential must have charge 2 and therefore is 
an odd polynomial, 

N N 

W{<^') = + - ^iifc*'*^*''- (2-35) 

3 

4=1 i,j,k=l 

Such an R-symmetry is clearly non-generic. In particular, the origin in these coordinates 
must be a special point for such a symmetry to hold, as an expansion of the superpotential at 

^By generic, we mean that a certain polynomial in the coefficients known as a resultant is non- vanishing 
-see ami 

^As a test, we note that the prototypical example of a landscape superpotential, the Gukov-Vafa-Witten 
superpotential for the complex structure moduli of a IIB orientifold model indeed scales this way - see 
1)2. 28(1 . Again, in a appropriate basis, the nonzero coefficients in the Gukov-Vafa-Witten superpotential can 
be ~ C(l), but, they are, indeed, very sparse. 

'^See |28I2HI8()I81I82I88I84I85| for discussions regarding the scale of supersymmetry breaking in landscape 
models. 
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any nearby point in field space certainly would not exhibit the same R-symmetry. That is, the 
R-symmetry is spontaneously broken by a vev for the In fact, given an arbitrary SUSY 
non-linear sigma model, there is no reason to believe that its superpotential will generically 
ever have a point in the target space where W has such a symmetry. Thus, the imposition 
of an R-symmetry means that we are restricting our consideration to very special SUSY 
non-linear sigma models expanded locally about a special point in their target spaces. This 
is the price that one must pay for a low-energy explanation for the scanning of the vacuum 
energy. For further discussion regarding this issue, see |36I37I38'] . Of course, R-symmetry is 
phenomenologically desirable for many other reasons as well (see the discussion in jlOj). 

The algebraic consequences of the R-symmetry and their geometric interpretations will be 
of use in our discussion of the statistics of the vacua of this model. Note that the conditions 
for unbroken SUSY are actually equations for the scalars 0* in the chiral superfields and 
the R-symmetry acts as a Z2 parity on these scalars, 0* — —0*. In particular, the SUSY 
vacua of this model are determined by quadratic equations, 

N 

d,W = A+J2 = (2.36) 

j,k=i 

which are invariant under the Z2 parity. As a result, these quadratics actually determine an 
ideal in the ring of Z2 invariants constructed from polynomials in the 0*, 

{diW) = {Ai + A,k(p'4>') e C[4>\ 0^]^^ (2.37) 

We can interpret the ring of invariants C[0^, . . . , 0^]^^ as the "polynomial functions" on the 
orbifold /Z2. Then, in analogy with the general case, we can consider the quotient ring, 

C[0\...,0^]^V(9,W^) ^ C2^"\ (2.38) 

which is a 2^~^-dimensional vector space over C generated by the images of the even monomi- 
als ■ ■ ■ <l'*^'' with ii < ■ ■ ■ < < N . Geometrically, this quotient ring can be interpreted 
as the polynomial functions on the variety cut out by the diW in the orbifold, which consists 
of the 2 images of the 2^ critical points of W in C^/Zs. 

2.5. Fermat Form of the Cubic 

We are interested in the properties of distributions of models with superpotentials respecting 
the R-symmetry ()2.35|) and quadratic Kahler potentials ()2.29|) . At first glance, one might 
expect to describe such a distribution of models as an arbitrary probability distribution on 
the space of all superpotential and Kahler couplings, f{Ai,Aijk,gij). However, we should 
not distinguish between models which differ by field redefinitions. Of course, arbitrary field 
redefinitions will not preserve the form of the Wess-Zumino models we are considering. 
However, it is easy to see that the field redefinitions which do are of the form 0* Gl^(f)"^, 
where G GL]\f{C). In particular, they act on the space of couplings as, 

Ai — >• AmG^, Aijk — >■ AmnpG^G^Gl., gij — >■ gmnG^Gj. (2.39) 
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If we posit some probability distribution on the space of couplings, this distribution gets 
averaged over G'Ljv(C) orbits to produce a probabihty distribution on the space of theories. 
That is awkward to deal with. If possible, it is much better to fix the redundancy by choosing 
a gauge for the GLn{C) field redefinitions. Since GLn{C) is dimensional, such a gauge 
condition should involve precisely iV^ independent, complex algebraic conditions. Now, in 
working with algebraic equations, it is often convenient to redefine our variables in order to 
make the coefficients of the terms of highest degree as simple as possible. If is a generic 
polynomial of degree d, we can use the GLj^{G) symmetry to enforce A^^ conditions on Aj^, 

Ai...i = 1, Aij...j = Ajij...j = ■■■ = Aj...ji = 0, for i ^ j. (2.40) 

We will refer to a polynomial satisfying these conditions as being in Fermat form, 

W = ^ {{(p^y + ■■■ + {(p^y) + (terms of degree < ci - 1 in each (p'). (2.41) 

(Jj 

Note that if W is odd (so it has a Z2 R-symmetry as above), then the Fermat form is 
somewhat stronger, as the additional terms in 1)2.411) actually have degree at most {d — 2) in 
any 0*. In particular, for the case of interest here, d = 3, with the R-symmetry above, the 
superpotential takes the form, 

TV N 

^i^') = E - - E b^,k<P'4>^<P^ (2.42) 

i=l i<j<k 

where a, has mass dimension 2 and is generically x 0{1) and the bij^ are symmetric, 
traceless (so bm = bijj = bjij = bjji = 0), dimensionless and generically \bijk\ ~ 0{N~^) 
with at most 0{N) of them 0{1). In going to Fermat form, it is important to note that we 
cannot also simultaneously set Qij = 6ij. Further, the genericity of W is important here — 
not all polynomials can be put into Fermat form. For example, the following model with 
spontaneous breaking of supersymmetry via the O'Raifeartaigh Mechanism, 

W = Z{X^ -a) + YX^ (2.43) 

cannot be put into this form. In fact, it is easy to show by analytic computations using 
eigenvalue methods (see Q that for = 3 and d = 3 no polynomial in Fermat form 
exhibits the O'Raifeartaigh Mechanism. 

We believe that this may be true much more generally. More precisely, we conjecture 
that all cubic, R-symmetric superpotentials of the Fermat form ()2.42p have at least a pair 
of supersymmetric vacua. While a proof of this fact is far beyond the scope of this work, we 
hope to return to this question in the future. 

2.6. Symmetries of the Fermat Form 

Just as is often the case in gauge fixing, it is important to note that our gauge condition 
(12.40)1 does not completely fix the GLisi{C) redundancy. For any generic Ai^, there is a finite 
subgroup H of GLAr(C) transformations which transforms the coefficients in such a way as 
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to keep Ai^ in the Fermat form ()2.41|) . Namely, /i™ G H are the solutions to A^^ algebraic 
equations of degree d in the A^^ complex components /i^, 

A^,...^^hr---hT' = l, Am,...m,hrhf'---hY'' = 0, for 2^ J. (2.44) 

If these were generic degree d equations, Bezout's theorem would tell us that they have d^'^ 
solutions. However, it turns out that this is not the case for generic Ai^, and d^^ is actually 
a strict upper bound on the order of H. We will show this by exhibiting a large subgroup 
of H which is independent of Ai^ and whose order does not divide d^ . First, note that 
permutations of the 0* certainly won't take us out of Fermat form, so we expect that the 
symmetric group on variables C H. Next, note that multiplying any 0* by a rf*'* root 
of unity also won't ruin the special form, so we get a subgroup of H isomorphic to (Z^)^ 
given by, 

}4 = 5iC', i = e'-^'l\ l<ni<d. (2.45) 

Thus, we see that x (^d)^, which has order N\ x d^ , must be a subgroup of H. Since A^! 
generally does not divide d^ , we see that, 

Nlxd^ < \H\ < d^' (2.46) 

In the simplest non-trivial case, A^ = 3 and d = 3, we numerically solved ()2.44|) using 
Mathematica for various values of the coefficients and found that H is generically of order 
648 = 2^ X 3^ = 3! X 3'^ X 4. In particular, there are indeed elements of H which depend on 
Ai^, which suggests that \H\ may generally be strictly larger than Sn x (Z^)^. This is not 
unexpected. Rather generally, the space of common solutions of these polynomial equations 
is an algebraic variety. Different values for the coefficients of the polynomials will nonetheless 
lead to isomorphic algebraic varieties. H is the modular group, and we will not attempt to 
give a full characterization of H beyond its "obvious" S'at x (Z^)^ subgroup. We will content 
ourselves with exploiting the constraints that stem from this latter subgroup. 

To summarize, after fixing gauge and imposing R-symmetry, each model is uniquely 
described up to a discrete symmetry H containing Sn x (Z^)^ by a cubic superpotential in 
Fermat form ()2.42j) and a quadratic Kahler potential ()2.29p . Thus, a distribution of field 
theory landscape sectors is a probability distribution on the space of couplings, 

/(oi, bijk, gij)daidbijkdgij (2.47) 

which is invariant under H D x (Zs)^, 

Gj ^ (imhi , bij}^ > b^nphi hjh'^, gij > gmnhi hj . (2.48) 

We assume that this distribution is fixed by high energy physics. Lacking any clear physical 
input which distinguishes among the different vacua of a given low-energy model, we simply 
treat these vacua democratically. In the following sections, we use algebraic methods to 
compute some statistical properties of various physical quantities averaged over the 2^ vacua 
of each model as a function of the couplings. Of course, these model averages can then be 
further averaged over the ensemble using the high energy distribution function. 
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3. Eigenvalue Methods for Solving Algebraic Equations 

The solution of simultaneous linear equations in variables is something we all learn to 
do from an early age using several different methods. For example, we may decide to do 
it via Gaussian elimination or via matrix methods - using determinants and Cramer's rule 
to invert the matrix of coefficients. A natural question to ask is if there exist analogues 
of these methods which could be used to solve simultaneous equations of higher degree. 
Indeed, this is the case. The method of Gaussian elimination has a vast generalization in 
the study of algebraic elimination theory and Grobner bases methods (see |2(ip27j ). These 
methods provide algorithms which allow one to systematically eliminate variables one by one 
in any given system of equations at the cost of increasing the degrees of the resulting system 
of equations in fewer variables and, with some further work, find (approximate) solutions. 
Unfortunately, these are not directly useful for us as we will be interested in the statistical 
properties (such as the average and variance of the superpotential and other observables) of 
such solutions as a function of the coefficients - i.e. in the properties of solutions of families 
of such equations. It turns out that the generalization of matrix methods to the case of 
simultaneous equations of higher degree does end up being quite useful for these purposes. 

3.1. Rings, Ideals and a Toy Model - Quadratic W 

We will introduce these methods by applying them to the simple case of a generic quadratic 
superpotential, 

W = Wo- Azi + ^CijZiZj = Wo-A-z + ^z-C -z. (3.1) 

Of course, the "statistics of vacua" in this case may seem to be, well, vacuous - the critical 
point of this polynomial is given by the unique solution of A^ linear equations, 

diW = -Ai + CijZj = O^C-z = A, (3.2) 

which is, of course, 

z = C-^-A, (3.3) 

where the genericity condition is that det C 7^ 0, or C G GLjv(C). However, if we are given a 
distribution of coefficients for W, we might be interested in the corresponding distribution of 
the values of W at the critical point. To attack such problems more generally, it is essential 
(though, in this case, akin to trying to kill an ant with a machine gun) to recast these 
questions in the algebraic language of polynomial rings and ideals. 

To begin with, let us quickly review what we will need about ideals, rings, and varieties. 
Roughly, a ring is a set R whose algebraic properties mimic many of those of the integers 
Z. That is, we may add, subtract, and multiply elements of R and obtain new elements in 
R, though the same may not be true for division. Recall that an ideal / contained in a ring 
R is defined by the property that if /, (7 G / and r,s ^ R then rf + sg ^ I. In particular, 
this implies that the quotient ring oi Rhj I, Q = R/I where we identify all elements of R 
differing by an element in /, is well-defined. Note that all the rings that we consider are 
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Noetherian rings, which means that all their ideals are finitely generated, so 



I={h,...Jra) = {j2 ^*/^ I ^« e ^} = {f^)■ (3-4) 

The rings we will be most interested are rings of polynomial functions in variables, R = 
C[zi, . . . , zn], and their quotients by ideals / generated by some polynomials fi{zi, . . . , z^) G 
R. Thus, I is just the set of all polynomials which vanish on the zero locus in of the 
simultaneous set of polynomial equations fi{zj) = 0, i = 1, . . . , m. This zero-locus is called 
the algebraic variety in corresponding to the ideal I. Thus, non-vanishing elements of 
the quotient ring Q = R/I are precisely the non-trivial residues of polynomial functions on 

restricted to the variety - they define a notion of "polynomial functions" on the variety. 

As we have already mentioned above, the polynomials diW which determine the critical 
point generate an ideal in the ring of polynomials in variables, 

/ = {diW) = {-A, + Cijz,) = {-A\C-z). (3.5) 

Clearly, any linear combination of the generators with complex coefficients is also in 
the ideal. In particular, if det C 7^ then C is invertible and we can consider the A^ elements 
of /, 2; — ■ A, obtained by multiplying the A^ generators diW by the matrix C~^. In 
fact, as C is invertible, these are also generators of the ideal /. Thus, in the quotient ring 
Q = C[zi, . . . , zj^]/!, we can use the N relations Zi — C^^ijAj = to systematically eliminate 
the Zi and rewrite any element of Q in terms of multiples of the identity - constants. So, we 
see that 

g = C[^i,...,^jv]// = C, (3.6) 

a one dimensional complex vector space, generated by constants. In particular, the residue of 
any polynomial f{zi, . . . , zn) G C[zi, . . . , z^] in Q is obtained by setting to zero all elements 
of the ideal /, that is, by substituting z = - A G into /. Thus, we see that the residue 
of f in Q is precisely its value at the critical point ofW, and that the variety corresponding 
to / is the single critical point z = ■ A G C^. Indeed, as the only functions on a point 
are constant functions - the value at the point - this is consistent with our intuition that the 
quotient ring Q = C[zi, . . . ,zn]/I should be interpreted as the "polynomial functions" on 
a point. Further, note that as Q is isomorphic as a ring to C, the product structure of the 
ring of polynomials is preserved in the quotient. That is, we can think of the residue of / in 
Q as an operator on other elements of Q which is multiplication by its value at the critical 
point. In particular, using the fact that dW, z — ■ A G /, we can compute the residue of 
W in Q by, 

W = Wo - ^A ■ z + ^z ■ dW Wo - ^A ■ z Wo - ^A ■ C-^ ■ A e Q. (3.7) 

Of course, we could have just plugged in our solution, but it turns out that this kind of 
computation generalizes to the case of higher degree in a way that is useful for computing 
statistical properties. 
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3.1.1. N on- Generic Quadratic Superpotentials 

It is also interesting to ask what happens if C is not generic, if detC = 0. Then, C has a 
non-trivial null space which is the space solutions of the homogeneous linear equations we 
get in the limit that we take Ai 0. We can describe this limit more precisely by lifting 
the equations to projective space, to CP^. That is, we add a new coordinate zq (which we 
may think of physically as a chiral superfield corresponding to a (complexified) scale factor) 
and then consider the homogeneous system of N linear equations C^Zj — AiZo = up to 
(complex) scaling in CP^. In the open set Uq = {zi ~ Xzi G C^+^-{0}\zo ^ 0} C CP^, the 
Ui = are good coordinates and the above equation reduces to the inhomogeneous equation 
with Zi Ui. Thus, if detC ^ 0, the unique solution of the inhomogeneous equation in 
Uo gives us the unique solution in CP^ as well. However, the advantage of projectivizing 
the problem is that we have added points corresponding to solutions even in the non-generic 
case. If det C = then we have at least one solution in CP^ which has zq = given by a 
one-dimensional subspace of the null space of C. As the — > oo as we take 2:9 0, we can 
think of such a solution as "a solution at infinity" from the perspective of Uq. Thus, we see 
that non-generic points in the space of couplings of the quadratic superpotential correspond 
to situations in which the supersymmetric vacuum has "run off to infinity" in field space. 
Now, since all of our field theory computations are approximations only valid locally in field 
space, these non-generic points should more properly be understood physically as points in 
the space of couplings in which those approximations break down. 

3.2. The General Case 

Now, let's consider the generalization of the above to the case that W is of degree d > 2, 

W{zu...,zn) = ^^n---^*.- (3.8) 

n=l,I„ 

The critical points of W are then given by the simultaneous solutions of degree d — 1 
polynomials, diW = 0, and, as we have already mentioned above, the polynomials diW 
generate an ideal in the ring of polynomials in variables, 

d 

I={dW) = { (3.9) 

n=l,I„_i 

Just as before, we may now try to simplify the generators by taking C-linear combinations 
of them. In particular, it would be nice to simplify the form of the terms of highest degree 
as we did in obtaining the Fermat form. So, define the matrix C by, 

Ai...i = Cii, Aij...j = Ajij...j = ■ ■ ■ = Aj...ji = Cij, for i ^ j. (3.10) 

Now, just as in the case of the quadratic superpotential, if det C 7^ ^ C G GLn{C), we 
can multiply the diW by the and obtain a simpler set of generators, 

d N 

I = { Y <^ij^^iin-i^n ■ ■ ■ ^iu-i) = (4'^ + (terms of degree < d - 2 in any zj)). 

n=l,I„_i j=l 

(3.11) 
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Of course, if W is in Fermat form, then Cij = 6ij and so the generators diW themselves 
take precisely this form.^ As we have argued earlier, the Fermat Form is convenient for 
other reasons as well, so we will assume henceforth that W is given in Fermat form. Just as 
before, the above form suggests that in the quotient ring Q = C[zi, . . . , zn]/I we should be 
able to use the generators in (jH.llj) to systematically eliminate all powers of Zi greater than 
{d — 2) and write any element of Q only in terms of monomials with powers of Zi less than or 
equal to (d— 2). This is in fact the case as long as TV is a generic polynomial of degree d (for a 
rigorous proof of this fact, see Chapter 3, Theorem 6.2 of [2^1 )• It is not difficult to show this 
explicitly in the case of interest d = 3, but the result is not particularly enlightening and will 
not be presented here. We only note here that the computation requires that several large 
matrices in the coefficients be invertible, which we believe is a computational manifestation 
of the assumption of genericity. For the following, we will proceed assuming that W is 
generic, and discuss what we mean by this more precisely at the end of the discussion. 

Thus, for W generic, the monomials zf^ ■ ■ ■ 2;^ with < di < d — 2 form a basis for Q 
as a vector space, and the relations will allow us to express the residue of any polynomial 
f{zi, . . . , zn) G C[zi, . . . , zn] in Q in terms of this basis. This is particularly simple in the 
case of the superpotential itself, as we have 

1 d.— l ^ d—l ^ 

^=d^ + ■ ■ ■ ^ Yl ■ ■ ■ e Q. (3.12) 

i n=0,I„ n=0,l„ 

Thus, we see that the quotient ring, 

C[z,,...,ZN]/{d,W)^C^''''^\ (3.13) 

is a (c? — 1) ^-dimensional vector space over C generated by the images of the monomials 
Zi^ ■ ■ ■ z'l^ with Q < di < d — 1. Further, as we mentioned earlier, this quotient ring should 
be understood as the "polynomial functions" on the variety cut out by the generators of 
the ideal, which in this case are the functions on a set of (rf — 1)^ points. In particular, 
in analogy with with the quadratic case, we expect that the residue of any polynomial 
f{zi, . . . , zm) € C[2;i, . . . , z^] in Q should be related to the values of / at the critical points 
of W . However, the fundamental difference between the case d > 2 and the simple case 
of a quadratic superpotential is that here, the residue of f in Q is a vector, and it is not 
immediately obvious how the components of that vector in any given basis are related to its 
values at the critical points. However, as Q is a ring, / also acts by multiplication on Q. 
Since Q is a vector space, and as multiplication by / is C-linear, / is also a linear operator 
on Q, which we can represent as a matrix f in our basis of monomials in an obvious way. In 
particular, we could consider / operator or matrix Zi in this sense. Clearly, the 

action of Zi on most basis vectors (for d^ < d — 2) is totally obvious, 

Zi : ■ ■ ■ zf ■■■z^" ^ ■ ■ ■ +^ ■ ■ ■ 4^, d,<d-2. (3.14) 



^Now, while it is obvious that a superpotential in Format form has dctC 7^ 0, one might wonder if this 
condition is in general necessary and sufficient for the existence of a Gijv(C) coordinate transformation 
taking any given superpotential to Fermat form. However, this is not at all clear, as C involves only nearly 
diagonal components of Ai^ while the transformation to Fermat form H2.39|l involves all the components of 
Ai^. We will not need to delve into this issue further for our purposes. 
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However, for di = d — 2 we must use the relations in the ideal to re-express the resulting 
monomial as a linear combination of the basis elements. Assuming that this is done for 
each Zi, we obtain N matrices Zi, each representing multiplication by one of the Zi. We can 
consider the characteristic equation for each of these matrices, 



which is a polynomial equation of degree {d — 1)^ in A. Generically, this equation will have 
{d — 1)^ distinct roots, corresponding to the {d — 1)^ eigenvalues of the matrix Zi. In partic- 
ular, the matrix Zi must be diagonal in the corresponding basis of eigenvectors. Since Q is a 
commutative ring, all matrices corresponding to multiplication by functions must mutually 
commute, and this must in particularly be true for the matrices Zi. Thus, all the matrices 
corresponding to multiphcation by functions / G C[zi, . . . , z^] can be simultaneously diago- 
nalized in the basis of eigenvectors of the Zi. In particular, in this eigenbasis, Q splits as a 
ring into the direct sum of {d — 1)^ rings isomorphic C, 



each of which one may naturally associate with the ring of functions on one of the critical 
points of W. Thus, generalizing our result from the toy model, the eigenvalue of f associated 
with each eigenvector corresponds to the value of / at the corresponding critical point. If 
we take f = Zi, then the eigenvalues of Zi are precisely the Zi coordinates of the critical 
points. Further, this means that the trace of f in any basis is the sum of the values of / at 
the critical points of W. In particular, we can do this trace explicitly in the monomial basis. 
Now, as powers of / correspond to powers of the corresponding matrix f, we can take their 
traces to compute the sums of powers of the values of / at the critical points as well. Thus, 
we see that by computing traces in the monomial basis, we can compute all holomorphic 
moments of / at the critical points of a generic superpotential W. In other words, we can 
effectively compute all (holomorphic) statistical properties of the values of any polynomial 
/ at the critical points of W. 

3.2.1. Resultants and Genericity 

There are two basic ways in which the generic situation of {d — 1)^ vacua can break down. 
Roots of the system of polynomials can coincide, or roots can run off to infinity. The former 
is something that already was implicit in the setup of ^0], where the chiral fields were 
assumed to be decoupled. When roots run together, the tunneling between the respective 
vacua is no longer suppressed, and one should not really count them as independent vacua. 
It is not difficult to modify the following considerations to deal with this situation, but we 
will not discuss this further here (see j2Sl for details). Having roots "run off to infinity" is 
a phenomenon that can occur only when interactions between the fields are turned on, and 
our analysis breaks down in this case. 

It turns out that both the collision and expulsion of roots can be captured by a single 
genericity condition.^ As it is only the latter non-genericity which is fatal, it is useful to have 

^It is a certain integral polynomial in the coefficients Ai^ of W known as the u-rcsultant. 




(3.15) 



Q = C[Zi, ZN]/{^^W) ^ C © ■ ■ ■ © C . 



(3.16) 



(d-i) 
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a precise criterion when it obtains. In the case of a quadratic superpotential, we saw that the 
critical point runs off to infinity precisely when detC = 0. Just as before, we can describe 
this more precisely by considering the homogenization of the polynomial equations dtW = 
and lifting them to to equations in CP^. It turns out (see jSHj Chapter 3) that there is an 
analogous integral polynomial in the coefficients of the monomials in W of highest degree Ai^ 
known as a multi-polynomial resultant which plays the same role as det C in degree > 2. The 
vanishing of this particular resultant indicates that some solutions have "run off to infinity" , 
or equivalently, the presence of solutions in the compliment of the open set Uq. We will refer 
to the vanishing locus of the resultant as the discriminant locus. While we will not describe 
this resultant in general, we will be able to compute it exactly in an example d = N = 3 
and see the advertised behavior explicitly in the next section. Finally, note that the cubic 
polynomial superpotentials we are interested in are approximations valid only within some 
small polydisk about the origin. Indeed, in the real situation, we don't literally have roots 
running off to infinity. Once they leave the polydisk we are considering, our approximation 
of truncating to the renormalizable terms in the Lagrangian breaks down. We simply no 
longer trust those solutions which have wandered too far from the origin. 



4. Holomorphic Moments and the Statistics of SUSY Vacua 

We will now restrict our consideration to the physically relevant case of a renormalizable, 
R-symmetric, cubic superpotential for TV chiral fields in Fermat form, 

N N 

W{zi) = ^ {^zf - ttiZi) - ^ bijkZiZjZk, (4.1) 

i=l i<j<k 

where we recall that bijk is symmetric with vanishing diagonals, i.e. bm = buj = biji = bjn = 
0. Dimensional analysis as well as the x (Zis)^ symmetry will be important tools for the 
analysis to follow. As the superpotential W has mass dimension 3, the Zi must have mass 
dimension 1 while the Cj have mass dimension 2 and the bijk are dimensionless. Further, 
C G (^s)^ acts as, 

Zi ZiCi, ai a,Cr\ hjk bijkC^Ci^Ck^, C = {0}, (4-2) 

where Q is in the z*^ Z3, under which Zi has charge +1 while the Oj and bijk have charge — 1. 
Note that the equations for the critical points of the cubic in Fermat form are, 

Zi - CLi-^ bijkZjZk = 0, (4.3) 
j<k 

and a basis for the quotient ring Q is given by the 2^ monomials, 

{l,Zi,ZiZj,...,zi---ZN}, (4.4) 

where at most one power of each Zi appears in each monomial. The residue of in Q 
written in this basis is. 
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Note that tr ( W) and all other holomorphic moments of odd powers of W or Zi vanish due to 
the R-symmetry. Explicitly, since the relations ()4.3p are even, any odd power of the operator 
Zi maps monomials of even degree to those of odd degree and vice- versa and therefore never 
has any non-zero diagonal components. 

We will be interested in determining the distribution of values of W among the 2^ vacua 
which appear for generic values of the couplings. These are the eigenvalues of W. By the 
R-symmetry, the values of W at the critical points occur in pairs, ±A, and are the roots of 
the characteristic equation, 

2N-1 

= det(|W - XI) = P{X') = {Xy-" - 2^-1 Yl (4.6) 

k=l 

We can express the Fk{a, b) in terms of the holomorphic moments of W by noting that. 



(2N-1 
- V ^ ( -"T" ) 
k=i 2fc V2 A y 



2k 



(4.7) 



where the last equality is only meaningful for positive powers of A and encodes the form of 
the Fk{a,b). Dimensional analysis and invariance under the Sn k (Za)^ symmetry can be 
used to restrict the form of the Fk{a, b). They must be homogeneous polynomials in the a^, 
of degree 3k, whose coefficients are rational functions of the bijk- Of particular importance to 
us is Fi, which is proportional to the "holomorphic variance" of W. Using the symmetries, 
we may parameterize this variance as, 

((|W2)) = 2-^tr(|W2) =Fi(a,6) = /(6)^a3+ ^ g^^\b)a,a,a, + Y^h^\b)a^a,, (4.8) 

i i<j<k i^j 

where the normalization factor converts the trace into the average value of ^W'^ among the 
2-'^ critical points. Now, as Fi must be invariant under ix (^3)^, we see that f{b) must 
also be invariant, while g^^^{b) is symmetric with vanishing diagonals and must have charge 
-fl under the i^'^, j*^, and k^^ Z3, and /i*-'(fe) has vanishing diagonal and must have charge —1 
under the z*^ and +1 under the j^^ Z3. While we will not be able to compute the coefficient 
functions /(&), g'^^^ip), and h^^ip) in closed form for arbitrary large A^, we can certainly do 
so for any small fixed TV. We will now turn to the simple example of the case = 3. 

4.1. Example I - W odd, = 3, iV = 3 

Consider the case of the generic, odd, cubic superpotential in three variables in Fermat form, 

^ = \ {4 + 4 + 4) - + Ci2Z2 + 032:3 + bZiZ2Z^) . (4.9) 



^'^Note that in general, the number of roots that run off to infinity as one approaches the discriminant 
locus is controlled by the order of the pole of F2N-1 (a, b). For a pole of order m, 2m roots run off to infinity. 
If the order of the pole is less than the maximal (2^~^), then some of the roots remain finite, even in the 
limit. 



19 



a2 - bziz^ = 0, 



The equations for its critical points are, 
zf — ai — bz2Zs = 0, z^ - 
A basis for the quotient ring Q is given by the 2'^ monomials, 

{1, Zi, Z2, Zs, Z1Z2, Z1Z3, Z2Z3, ZiZ2Zs}- 

In this basis, the residue of in Q is, 

2 

W = -- (fli^i + 022:2 + as^s) ■ 
It is an easy exercise to determine the form of the matrix W, 



as - bziZ2 = 0. 
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fact that we can just as easily compute, 
((fW^)) = 2-hT (|W^) = 



if + a| + a|) (1 - jb^) + |aia2a362 
1-63 



ig is the 
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as well as all higher moments of W. In particular, we can explicitly read off the coefficient 
functions of ()4.8|) for the = 3 case, 

14-6^ 

fib) 



gib) 

Kb) 



Al-b^' 

9b' 



(4.15) 
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Further, if we wish to explicitly compute the values of W at the critical points, we can use 
these moments to compute the characteristic polynomial for W and numerically solve it. 
Indeed, this is easy to do using Mathematica, and one can verify that numerical solution of 
the above equations through other means agree with the eigenvalue methods. Further, using 
similar techniques, we can also compute the moments of any polynomial over these vacua. 

Finally, we note that ()4.15|) are singular in the limit that b approaches a third root of 
unity. It is easy to check numerically that six of the eight critical points run off to infinity in 
this limit, and that these three points indeed comprise the discriminant locus of the = 3 
case. Thus, we see explicitly in this case that the multi-polynomial resultant we discussed 
earlier must be proportional to (1 — b^) for = 3. 
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4.2. Example II - The Fermat Cubic for Large and small bijk 

Since the dimension of the vector space Q grows exponentially with the number of variables, 
the matrix methods introduced above become quickly intractable, even numerically, for 
> 10. In particular, explicit computations require that we can compute holomorphic 
moments like 

((fW^)) = 2-^tr (fW^) = ^a,a,(2-^tr (z^Zj)) = J] a.a,((ziZj)), (4.16) 

which in turn require that we can compute the moments of products of the Zi's. This could 
be done if we could easily compute the matrix elements of Zi in this basis. However, as 
we mentioned earlier, this is very computationally intensive, and involves the inversion of 
exponentially large matrices. 

Thus, it is useful to understand if there exists a limit in which the above methods become 
more tractable at large A^, at least perturbatively. Now, note that the decoupled limit with 
bijk = considered in [Hlj is certainly such a simple case, so the bijk seem to be natural 
candidates for small parameters. Indeed, radiative stability of the Kahler potential requires 
that generically \Aijk\ ~ 0{N~^) while at most 0{N) of them may be 0{1). That is, the 
bijkS are generically small. There are, of course, a large number of them, so this large-iV 
suppression of the magnitudes of the bijk is compensated by the sum. But say that we make 
the further assumption that perturbation theory is good — that is, that the magnitude of 
the bijk is further suppressed by some small parameter 6. Then it makes sense to expand 
our expressions for ((/(z))) as a (formal) power series in the bijk- 

4-2.1. Computations of Holomorphic Moments for Small b^k 

To begin with, we can use dimensional analysis as well as the Sn x (^s)^ symmetry to 
constrain the form of successive terms in the power series expansions for ((/(z))). Let us 
consider ((ziZj)) as an example. First note that if we have i = j, 



(fljl + ^ feijfcZjZk)) = + ^ 6ijfe((zjZk)), (4.17) 

j<k j<k 



SO we only need focus on the case i ^ j- Using the symmetries and dimensional analysis, it 
is easy to check that the first few terms in the (formal) power series expansion of ((zjZj)) for 
i j must take the form, 

((ziZj)) = af ^ ^ + ^ bikibjki af^ {bjkiai + hkiaj) + ^ [hjkai + hjiak) + 0{b'^), 

k k<l 

(4.18) 

where the ^ e Q are rational numbers which we will see are independent of TV for any 
fixed a^f^ and sufficiently large. The a^f^ can be found by computing the diagonal matrix 
elements of the operator ZiZj in the monomial basis, summing them, and dividing by 2^. 
Of course, the only way one obtains a diagonal element is if the relations ()4.3p are used, 
so such diagonal elements only arise from the action of ZiZj on basis elements containing Zi 
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or Zj. Each use of the relation results either in the introduction of one factor of or 

so the a-'^'' can be computed using the relations precisely (ri + 1) times. For example, the 

contribution of the monomial Zi to af'^ is found by, 



= h ^ hi.ji{ajZi+ ^ hjmnZmZnZl+ ^ bjimZ\Zm ' Zl^ 

m<nj^{l,j) rnj^(j,l) -j^g^ 

^ ^ 5Z 5Z ^ijl^jlrnaiZm + 

= ■ ■ ■ + + ■ ■ ■ • 

k 

In fact, we can greatly simplify this computation by noting that the Sn symmetry implies 
that we only need to compute how many times any single term of a sum over dummy indices 

(2) 

appears in the diagonal matrix elements. For example, to find , we only need to ask how 
many times the single term b'^jj^ak with i,j, and k all fixed appears on the diagonal of ZiZj. 
This is very easy to do. First of all, this term only involves ak and one particular coefficient 
bijk- As these two coefficients arise in the relations ()4.H|1 only in terms involving Zi, zj, and 
Zk, we can completely ignore all the other variables, which just go along for the ride.^^ Thus, 
to calculate a^i \ we only need to compute the coefficient of the term b^-^^ak in the trace of 
ZiZj restricted to the 2^ monomials of the form z'l^z'^^z'^^ , a, = 0, 1 in our basis ()4.4|1 . ignoring 
all terms involving any other variables (and divide the result by 2'^). That is, we have shown 
that this coefficient can be computed by considering just the case = 3. For instance, we 
can expand ()4.14|) in a power series in b and find a{ 



(2) _ 3 



4' 



We can more easily compute this using a diagrammatic technique. To compute the 



contribution of ZiZj ■ to af^, first draw a labeled external line for Zi and zj as 

well as all the 2;„'s present in the basis element. Each bijk corresponds to a vertex which 
takes in two identically labeled lines corresponding to one of its indices and outputs two 
lines corresponding to the other two indices, while at is an external sink for a pair of Zk 
lines. Construct all possible graphs with two bijk vertices, and one ak external sink from the 
given lines, noting that graphs which differ by a choice of which lines of the same index are 
contracted at a vertex or sink are equivalent. Count the number of distinct graphs you drew 
and divide by 2^ - this is the contribution. It turns out that there is precisely one possible 
graph for each basis element with either Zi or Zj in it, so we again find a\ 



(2) _ 3 



4' 

The above method, of course, generalizes in the obvious way to the computation of any 
af''^ by restricting ourselves to just the relevant = n + 1 variables. For example, to 
compute the coefficient a^f^ of bikib^ki^i^ consider A^ = 4 and graphs with one biki vertex 
and two bjki vertices as well as a single sink Oj. Note that it is useful to start with monomials 
of lowest total degree and proceed to those of higher degree as any graph associated with 
a given basis monomial also contributes to all monomials it divides. Thus, without much 

^^Formally, this is equivalent to working in Q modulo the ideal (z;), I ^ i,j, k. 
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trouble, one can compute the terms cubic in b. 



'ijki 



^i^j)) = 4 E ^'ijk"'^ + X] ^ikibjki [{hjkiai + hkiaj) + 2{hijkai + hjiak)] + 0{h'^) , (4.20) 



k<l 



SO ^ = |, a'^^ = |. Using this result as well as (|4.17p we can compute, 

((fW^)) = ($^a^((zi^)) + 2$^a,a,((z:Zj)) J = + ( 2a,a, + 6,,, | {{z,z^j 

\ i i<j / i i<j \ I 

1 + ! 5Z ) 5Z «^ + Yl 1 ( ^Ifc + 2 5Z bijibjkikki 1 a^a.-afc 

V i<fc</ / i i<j<k \ I / 



(4.21) 



Thus, we can read off the coefficients functions of ()4.8j) up to terms of order b^, 

j<k<l 

9''\b) = I (^bl, + 2 J2 b^,lb,kAk}j +..., (4.22) 
h^^{b) = IJ2^.kib%, + .... 

k<l 

As noted, each additional power of b brings an additional sum over a dummy index, which 
has 0{N) terms. If \bijk\ ~ N~^5, then the term of order b^ goes like A^^'^, so this is a 
systematic expansion in powers of 5. 

We can continue and compute higher holomorphic moments in an expansion in powers 
of bijk - for instance, the following holomorphic moments will be useful (with i j ^ k ^ I), 

i^,'))-W)f = Y^lka,ak + 0{b^) 



j<k 



(Zi^Zj))-((Zi^))((ZiZj 



■J2blk(^^^k + 0{b') 

k 

Zi^Zj^)) - ((Zi2))((zj2)) = Y^hkib.kiakai + (9(6' 



k<l 



(zi^Zj^)) - ((ziZj))^ = a-iaj + ^bikibjkiakai + 0{b^) 



(4.23) 



k<i 



((zi^ZjZk)) - ((zi2))((zjZk)) = hjkajak + 0{b^) 

((zi^ZjZk)) - ((ziZj))((z;Zk)) = bijkttjak + fj^ ^jfc;«i«/ + ^(b^) 

I 

JziZjZkZ,))-((ziZj))((zkZi)) = 0(6=^). 
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Therefore, we see that it is possible to effectively compute holomorphic moments and any 
desired holomorphic statistical property in the limit of large as a systematic expansion in 
small bijk- 

For any fixed N, we can, with more effort, crank out some explicit formula, capturing 
the dependence on the discriminant locus. For = 3, we see explicitly, that each of the 
moments in ()4.23|) has a double pole at the discriminant locus. 



b\8{2 + b^)a2a3 + Sb^aj) 
16(1-63)2 

3 2\\// \\ b^{{lQ-b^)aias + l2bal) 

Zi Z2 - Zi Z1Z2 = 



16(1-63)2 

/ 2 2\\ // 2\\// 2\\ b^{il6-b^)aia2 + l2bal) 

Zi Z2 )) - Zi )) Z2 )) = 



(zi^Z2^)) - ((Z1Z2' 



(zi2z2Z3)) - ((zi2))((z2Z3^ 



16(1 - 63)2 

2 _ 8(2 + b^)aia2 + 36^ai 
^ ~ 16(1-63)2 

6(8(2 + 63)0203 + 36^0^) 

16(1-63)2 
6((16 - 63)a2a3 + 126a?) 



(4.24) 



((zi2z2Z3)) - ((Z1Z2))((Z1Z3); ^^^^ _ ^3^^ 

5. Holomorphic Moments and Scanning 

Now, let us consider the relationship between the holomorphic moments we have computed 
and the real-valued statistical properties of the vacua. W = Wi + iW2 is complex. Averaged 
over the 2^ vacua, the mean value is, of course, zero. 

m)) = {{W2)) = 

The variances involve the sums of squares of the the eigenvalues, and so are encoded in 
|2Fi = ((W^)). The variances we are interested in are SW^, 6W2 and the covariance, 
(lWiW2))- Two linear combinations of these are "holomorphic" in the eigenvalues, 

6W^ - 5Wl = Re ((W^)) (5.1a) 
{{W1W2)) =lm{{W'~)) (5.1b) 

The remaining linear combination requires more detailed knowledge, though it is easy to 
find a lower bound from the Cauchy-Schwarz inequality, 

6W^ + 6W^ > \{{W^))\ (5.1c) 

In the following, we will turn to the implications of these formula when one averages over 
some ensemble of couplings. A lower bound, like ()5.1b ) will be quite sufficient to prove that 
a given coupling "scans" . But one requires something of an upper bound in order to prove 
that a coupling doesn't scan. These methods generalize straightforwardly to other coupling 
which depend holomorphically on the 0j (superpotential couplings for the Standard model 
fields, holomorphic gauge couplings, etc.). 
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5.1. Ensembles of Theories and the Scanning of A 



Let us return to ()4.8|1 . 



E 9'''ib)a,aja, + J2h''ib)a^a,\, (5.2) 

i i<j<k ij^j 

which encodes the variance of W among the 2^ vacua for fixed values of the couplings. What 
we would like to do now is discuss what happens if we have some ensemble of theories. For 
example, rather than studying F-theory vacua for fixed fiux, we might wish to study the 
ensemble of all possible fiuxes. 

No one, currently, has a compelling proposal for what probability measure to choose 
for that ensemble. Neither do we. What we hope to do in this section is to explore the 
implications one can extract from such a choice, given the results of the previous sections. 

The authors of ^U] worked with bijk = 0, and assumed that the were N independent 
random variables, chosen from some common distribution. In their limit, the second two 
terms of ()5.2|) vanish, and we have 

JV 

1=1 

Since there are terms, and each is an independent random variable, the sum grows like 
0{N). At first glance, for nonzero bijk, the situation appears to change dramatically. The 
second term in ()5.2j) contains 0{N^), and the third term contains 0{N^) independent random 
variables. At least, that would be the case if the bijk coefficients were all generically of 0{1). 
However, as we have argued, radiative stability requires that the bijk be generically small. 
So, despite appearances, each of the terms in ()4.8j) is 0{N). The cosmological constant scans 
(and, in fact, the variance grows like 0{N)) in this more general context, just as it did in 
the model of (TUj where bijk = 0. 



5.2. Other Couplings 

Other holomorphic couplings which depend on the (pi can be treated similarly. Consider 
some such coupling, c{(j)), which might be a holomorphic gauge coupling, or a coupling in 
the superpotential for the Standard Model fields. It is reasonable to assume that c(0) has 
definite parity under the Z4 R-symmetry that constrained the form of W{(j)). The crucial 
distinction in ^U] is between those cases where c{(j)) is odd and hence {{c{4>))) = (like the 
superpotential, W) versus those for which c(0) is even. The former couplings "scan," whereas 
the latter do not: the standard deviation of c(0) among the 2^ vacua is much smaller that 
its mean value. 

As with the superpotential, we assume that c{(j)) has a Taylor expansion, convergent in 
a polydisk of radius Mc- For simplicity, we will take it to be dimensionless; the case of the 
/i parameter in the Standard Model is an easy generalization. So we have, 

Cev(0) = Co + ^ CijCpiCpj/Ml + ... 

(5-3) 

Codd(0) = 2^ CiCpi/M^ + ... 
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depending on the parity of c(0). As in our previous discussion, radiative stability constrains 
the form of the coefficients in this expansion for large-A^. The Cijk,,. must be generically 
small. Specifically, we have constraints of the form 



fg^cifc^, ~ 0{N) 



det ii{g- 



3] 



0{l) 



(5.4) 



and so forth for the higher coefficients. Thus, the mean value of the (2/c)*^ term in the series 
()5.3j) (the mean value of the odd terms vanish), rather than going like N'^^ {Mr / Mc)^'' actually 
goes like N^Mr/Mc)"^^. So, we have a systematic expansion in powers of = {Mr/Mc)'^. For 
sufficiently small e, the leading, term can compensate for the overall factor of A^, while 
the subleading terms are negligible. More explicitly, if we let Oj = M^Oj we have. 



Codd) 







/WE 



(5.5a) 
(5.5b) 



where the 0{e^) terms represent the contributions of quartic and higher terms in c(0). Note 
that the rational functions of b: f{b), g'^^^{b) and K'^{b) are the same ones from ()4.8p that 
appeared in the computation of the variances of the superpotential and were computed up 
to terms of 0{b^) in (j4.22|) . It is easy to check that, as announced, the leading term behaves 
as Ne^, with the subleading terms suppressed by higher powers of e^. 
The variance is calculated similarly. In the odd case. 



-^odd/ 



+ 0(e^). 



(5.6) 



Thus, given a distribution for the couplings we can definitively show that this coupling scans 
if this variance is non- vanishing. 

In the even case, let c = c — ((c)). The variance is then given by. 



i,j,k,l 



(ziZjZkZ,))-((ziZj))((zkZ,))] + 0(e^ 



(5.7) 



As in §4.2^ we can easily compute ()5.7|1 in an expansion in powers of bijk using ()4.2H|1 . For 
finite A^, we can do better - for A^ = 3, we have ()4.24|) . As with the superpotential, we can 
write c = ci + ic2, and we have 



5cl - Scl = Re ((c') 
{{c,c,)) = Im ((c^) 
5cl + 6cl>\{{c'))\. 



(5.8a) 
(5.8b) 
(5.8c) 



The right-hand-side of ()5.5b ) goes like A^e^, whereas the right-hand-side of ()5.7p goes like 
A^e^. We would like to conclude that the standard deviations, Sci/ci and SC2/C2, behave as 
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Unfortunately, ()5.8b ) is only a lower bound, so more detailed information is necessary to 
really show that this coupling does not scan. Indeed, this is a potentially serious drawback 
to the whole notion of a "friendly landscape." If the lower bound in ()5.8b ) drastically under- 
estimates the true variance, then these couplings will vary appreciably over this ensemble of 
vacua. And anthropic arguments, based on holding them fixed while varying other coupling, 
like the cosmological constant, are incorrect. 

6. Generalizations 

The "space" of vacua discussed here is the complex affine algebraic variety, C[zi, . . . , zi\f]/ (diW) 
As such, it was amenable to the techniques of complex algebraic geometry. If we were study- 
ing real scalar fields, 0j, with potential, V{(f)), we would be faced with a problem in real 
algebraic geometry. This problem is harder, both because real algebraic geometry is harder, 
and less well-developed than the complex case and because the characterization of the de- 
sired space of vacua is more subtle. We are not interested in ]R[xi, . . . ,xn]/ {diV). We are 
only interested in minima of V, as opposed to all critical points. At large A^, "most" critical 
points of V are actually saddle points, and it's algebraically a little awkward to pick out just 
the minima. 

A more interesting generalization is the case in which some of our complex chiral multi- 
plets are charged under a U{1)'' gauge symmetry. Physically, this leaves open the possibility 
of supersymmetry-breaking, in this "moduli" sector, via the inclusion of Fayet-Iliopoulos 
terms. Mathematically, this gauging moves us from the realm of complex affine algebraic 
algebraic geometry to that of toric geometry. The toric version of our problem is nearly 
as well-developed mathematically as the affine case we have discussed. It would be very 
interesting to generalize our considerations to that case. 

Finally, we have deliberately eschewed discussion of the microphysics that determines the 
(ensemble of) couplings in our low-energy effective Lagrangian. Douglas and collaborators 
|39l4U|41l42l43l28l44|45l46l47j . for instance, have pursued the idea of treating all possible 
choices of fluxes in an F-theory compactiflcation "democratically," assigning equal weight 
to the low-energy theory that arises from each choice of flux. Further developments along 
this vein appear in ^48,.36,:49r5()f51f52p53j . It has been argued |22] that in the context of type 
IIA flux compactiflcations this may not be a reasonable choice, as there are instances where 
the number of possible choices of flux is inflnite. Independent of this more subtle question, 
we believe that our methods will prove useful in analyzing the vacuum structure of any give 
theory whenever one has a large number, A^, of light chiral multiplets. 
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